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ABSTRACT 

We have used high resolution (~2.3") observations of the local (D^ = 46 
Mpc) luminous infrared galaxy Arp 299 to map out the physical properties of 
the molecular gas which provides the fuel for its extreme star formation activ- 
ity. The 12 C0 J=3-2, 12 C0 J=2-l and 13 C0 J=2-l lines were observed with 
the Submillimeter Array and the short spacings of the 12 C0 J=2-l and J=3-2 
observations have been recovered using James Clerk Maxwell Telescope single 
dish observations. We use the radiative transfer code RADEX to estimate the 
physical properties (density, column density and temperature) of the different 
regions in this system. The RADEX solutions of the two galaxy nuclei, IC 694 
and NGC 3690, are consistent with a wide range of gas components, from warm 
moderately dense gas with T kin > 30 K and n(H 2 ) ~ 0.3 - 3 x 10 3 cm -3 to cold 
dense gas with T^ in ~ 10-30 K and n(H 2 ) > 3 x 10 3 cm" 3 . The overlap region 
is shown to have a better constrained solution with T^ in ~ 10-50 K and n(H 2 ) 
~ 1-30 x 10 3 cm -3 . We estimate the gas masses and star formation rates of 
each region in order to derive molecular gas depletion times. The depletion times 
of all regions (20-60 Myr) are found to be about 2 orders of magnitude lower 
than those of normal spiral galaxies. This rapid depletion time can probably be 
explained by a high fraction of dense gas on kiloparsec scales in Arp 299. We 
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estimate the CO-to-H 2 factor, a co to be 0.4 ± 0.3 (3 x 10" 4 / x C o) M Q (K km 
s -1 pc 2 ) -1 for the overlap region. This value agrees well with values determined 
previously for more advanced merger systems. 

Subject headings: galaxy: individual (Arp 299, NGC 3690, IC 694), infrared: 
galaxy 



1. Introduction 



Ultra/Luminous infrared galaxies (U/LIRGs) are systems in the local universe that 
exhibit extreme star formation. They emit a large portion of their total luminosity at far- 
infrared (FIR) wavelengths (LIRGs: L FIR ~ lO 11 " 12 L ; ULIRGs: L Fm > 10 12 L ; Sanders 
& Mirabel 1996). Early optical studies have shown that a large fractio n of U/LIRGs show 
morphologies that resemble systems that are interacting or merging (lArmus et al.l 119871 ; 
Sanders et al.lll988l ). High resolution optical and near-infrared imagi ng has revealed tha t 
a large fraction of U/LIRGs have nuclear separations of 0.3 - 48 kpc (IMurphy et al.lll996l ). 
The merger process likely triggers the extreme star formation in these galaxies (Sanders et 
al. 1988) but the detailed process is still not well understood. 

U/LIRGs have be en observed to conta in a large amount of molecular gas with M(H 2 ) 
up to about 10 10 M (jSanders et al.lll986[ ). The molecular gas is seen to b e concentrated 
near the nuclear regions within a radius of 0.5 kpc (jSanders fc Mirabellll996l ). Subarcsecond 
resolution CO obser vations have revealed ro tating disks of molecular gas that has been 
driven to the nuclei ( IDownes fc Solomonlll998l ). Theoretical models have shown that during 
the collision, the gas loses angular momentum due to dynamical f riction which causes the ga s 
to decouple from the stars and flow inwards towards the nuclei (IBarnes fc Hernquistlll99ll ). 
The concentrated molecular gas provides fuel for star formation and/or for an active galactic 
nucleus (AGN). 



Menendez-Delmestre et al.l (120091 ) have shown that submillimeter galaxies (SMGs) have 
similar radiation environments to local starbursts. The distribution of the optical depth 
determined from the silicate absorption feature at 9.7 fim was found to be lower than for 
local U/LIRGs, which suggests that the mid- infrared regions of SMGs have a lower dust 
obscuration. The difference in extinction between U/LIRGs and SMGs likely arises from 
differences in geometry. Using 12 CO J=3-2 observations, Iono et al. (2009; Paper II) have 
shown that the molecula r gas disks in local U/LIRG S are m ore compact (0.3 - 3.1 kpc) than 
the SMGs (3 - 16 kpc). iMenendez-Delmestre et al.l (120091 ) also found that the majority of 
SMGs are dominated by intense star formation just like local U/LIRGs. The differences 
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suggest that ULIRGs are not exact local analogs of high redshift SMGs; however, there are 
enough similarities that local ULIRGs are often used as templates for SMGs. 

Arp 299 (NGC 3690 + IC 694, Mrk 171, VV 118, IRAS 11257+5850) is one of the 
nearest (D^ = 46 Mpc) examples of a merger system. It has a far-infrared luminosity of 
Lfir = 5.5 x 10 11 L Q (Sanders et al., 2003), classifying it as a LIRG. In the optical, the two 
nuclei are still dis tinguishable. High resol ution observations show the nuclei are separated 
by 22" (4.9 kpc ; ISargent fc Scovillel Il99lh . Gehrz et al. (1983) showed that there is also 
a third region o f emis sion thought to be associated with the overlap of the two galaxies. 
Hibbard &: Yunl (Il999l ) found two 180 kpc HI tails as well as a faint optical tail displaced 
from the HI tails. They infer that the merger started 750 Myr ago and the two nuclei will 
completely merge in roughly 60 Myr. Evolutionary starburst models have shown that Arp 
299 has been going thr ough a recent episode of in teraction- induced star formation over a 
span of about 15 Myr (jAlonso-Herrero et al.ll2000l ). Interferometric observations of 12 CO 
show strong emission from the three major regions found in infrared images (Sargent & 
Scoville 1991; Aalto et al., 1997; Casoli et al., 1999) indicating the presence of large amounts 
of molecular gas. 

Studying the physical conditions and distribution of the warm, dense gas associated 
with star formation will help in understanding the processes and timescales controlling the 
star formation in galaxy mergers. In this paper, we analyze Submillimeter Array (SMA) 
interferometric observations of Arp 299 in the 12 CO J=2-l, 13 CO J=2-l and 12 CO J=3-2 
lines first published in Wilson et al. (2008; hereafter Paper I). In Section 2 we present 
the observations and the data reduction process including the short spacings correction. In 
§ 3, we present and discuss our radiative transfer modelling results as well the CO-to-H 2 
conversion factor, a co for the extranuclear region. In § 4, we discuss the star formation rates 
and depletion times for each region of Arp 299. 



2. Observations and Data Reduction 

2.1. Submillimeter Array 

We have used the SMA (Ho et al. 2004) observations centered on 12 CO J=2-l and 
12 CO J=3-2 for Arp 299 first published in Paper I. The correlator was configured to have a 
spectral resolution of 0.8125 MHz corresponding to roughly 1.1 km s _1 for 12 CO J=2-l and 
0.7 km s _1 for 12 CO J=3-2. The 12 CO J=2-l transition was centered in the upper sideband 
and 13 CO J=2-l was observed simultaneously in the lower sideband 10 GHz away. Since Arp 
299 has sufficiently extended emission, a small mosaic of two pointings was used. Further 
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details of the observations are given in Paper I. 

We start with the calibrated uv data from Paper I. These data were converted to FITS 
files and then exported into CASA format for further processing and imaging. The data were 
flagged to remove the first six and last six channels of each of the 24 chunks of the correlator. 
The data were also flagged to remove high amplitude value data points (removing less than 
1% of the data). The uv data sets were continuum subtracted using line-free channels. Clean 
data cubes with velocity resolutions of 20 km s _1 and 50 km s _1 were created using a robust 
weighting scheme in order to achieve a good balance between resolution and sensitivity. 
Each data set was cleaned down to 2 times the rms noise in each velocity channel (Table 
1). The beam was sufficiently clean and the emission is sufficiently compact that cleaning 
the entire inner quarter region produces good maps. However, for 12 CO J=2-l and J=3-2, 
we interactively placed clean boxes around emission regions in each channel. Integrated 
intensity maps were created with a 2a cutoff using velocity channels with emission from Arp 
299 (Figured]). Finally, the integrated intensity maps were corrected for the primary beam. 



2.2. Owens Valley Radio Observatory 

We have obtained a clean data cube of 12 CO J=l-0 observed using the Owens Valley 
Radio Observatory (OVRO) from S. Aalto first published in Aalto et al. (1997). Details of 
the observations and reduction are given in Aalto et al. (1997). An integrated intensity map 
was created with a 2a cutoff using velocity channels with emission from Arp 299 (see Figure 
1). The integrated intensity map was corrected for the primary beam. Adding this line to 
our analysis will help to better constrain the density and temperature of Arp 299. 



2.3. James Clerk Maxwell Telescope 

Arp 299 was observed with the James Clerk Maxwell Telescope (JCMT) on 2007 March 
31 for 12 CO J=3-2 (Program: M07AC11; PI: CD. Wilson) and 2011 March 09 for 12 CO 
J=2-l (Program: M11AC06; PI: K. Sliwa). The JCMT was used to make a 2' map with 
6" sampling using the array receiver HARP-B in the 12 CO J=3-2 line with a 15" beam and 
to make a 7" sampled 1.5' x 1' map with the RxA receiver centered on the 12 CO J=2-l 
line with a 22" beam. Further details of the 12 CO J=3-2 data are given in Paper I. The 
raw data were ma de into a cube tha t spanned from -400 km s _1 to 400 km s _1 using the 



Starlink software (ICurrie et al.ll2008l ) written for JCMT data. The cube was then baseline 



corrected using line free channels and a first order baseline. Several of the outer bad pixels 
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were trimmed off and the cube was binned to a velocity resolution of 20 km s l . 

We assume a main beam efficiency (r] mh ) of 0.6 for 12 CO J=3-2 and 0.69 for 12 CO J=2-l 
in order to convert the intrinsic units of antenna temperature (T£) to main beam temperature 
(Tinb = Tl/rj m h) and then to Jy km s _1 which gives us a scaling factor of 26.9 Jy K(T mb ) _1 
and 22.9 Jy K(T mb )- 1 for 12 CO J=3-2 and 12 CO J=2-l, respectively. We thus assume that 
emission from Arp 299 fills the beam, unlike Paper I where Arp 299 was treated like a point 
source. 



2.4. Short Spacing Correction 

We use the JCMT maps to recover the short spacings of the SMA maps. We combine 
the two maps in the image plane using the feather task in CASA. Both maps were binned 
to the same velocity resolution so that each velocity channel matches in each map. The 
feathering technique takes both maps, regrids the low resolution image to match the high 
resolution image, Fourier transforms them to the gridded uv-pl&ne, sums up the gridded 
visibilities and then Fourier transforms the feathered map back into the image plane (Figure 
2). Comparing the SMA only maps with the JCMT maps, the 12 CO J=2-l map was missing 
about 20% of the total flux and the 12 CO J=3-2 map was missing about 70% of the total 
flux (Table 2). Note Paper I estimated only 47% missing flux in 12 CO J=3-2 because they 
assumed Arp 299 is point-like in the K to Jy conversion of the single dish data. Casoli et 
al. (1999) have compared the OVRO 12 CO J=l-0 map of Aalto et al. (1997) with an IRAM 
30m single dish map and found that about 56% of the flux was missing. The missing flux 
in all the maps is an issue but we try to avoid the problem by using the peak integrated 
intensities for our analysis (see § 3.1). 

2.5. Line Ratio Maps 

For our analysis, we degraded the resolution of all the interferometric maps to 3.6" x 
2.4" in order to probe the molecular gas on similar physical scales. The elongated beam of 
the CO J=2-l maps limited the best possible degraded resolution to 3.6" x 2.4". We applied 
a Gaussian taper weighting to the SMA uv data to degrade the resolution. The resolution 
of the OVRO 12 CO J=l-0 map was degraded by applying a Gaussian kernel to the data 
cube. New integrated intensity maps were made from these resolution matched cubes. We 
converted the intrinsic units of the CO integrated intensity maps from Jy/beam km s _1 to 
K km s^ 1 . We then created 12 CO jE§Ef , 12 CO jEjz^ and J=2-l line ratio maps by 
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dividing the appropriate maps with a cutoff at 2a for each map (Figure [3]). Note the line 
ratio values near the edges of each region are noiser because the uncertainty by definition 
approaches ±50%. 



3. Radiative Transfer Modelling 
3.1. Physical Conditions in Arp 299 



We adopt the notation of iGehrz et al.l ( 119831 ) to denote each region of Arp 299. The 



two nuclei are denoted as region A for IC 694 and region B for NGC 3690. We denote the 
"disk" of region A as region A2. The strong emission to the north of region B, denoted as 
region C, can be broken up into two sub regions; the strongest peak to the west is denoted 
as CI and the peak to the east is denoted as C2 (Figure 2). Note that our sub-region C2 is 
the region C2+C3 in Aalto et al. (1997). 

The line width and size of the CO emission of each region have been measured using 
the 12 CO J=3-2 map which has the best angular resolution. The velocity line width FWHM 
of each region was fitted with a Gaussian and the results are presented in Table 3. The 
diameter of each region was obtained by fitting to a two dimensional Gaussian. For region 
C, each sub-region, CI and C2, was fitted to a two dimensional Gaussian because the entire 
region cannot be fitted well by a single Gaussian. The source sizes listed in Table 3 are 
deconvolved from the intrinsic 12 CO J=3-2 beam size. 

Using the line ratio maps, we obtain the line ratio values at the peak intensity pixel of 
each region found in the 12 CO J=3-2 map. The line ratios and la values are presented in 
Table 4. The la values are calculated assuming a 20% calibration uncertainty for each map 
for the 12 CO jZ^Zi an d 12 CO jZ^Zq li ne ratios and a 10% calibration uncertainty for each 
map for the jr^j J=2-l line ratio since these maps were observed simultaneously. We see 
that including short spacings with the JCMT maps does not change the peak ratio values 
significantly except for region A2. We assume that the missing flux for 12 CO J=l-0 and 
13 CO J=2-l at the peak intensity value is insignificant for regions A, B, CI and C2. The 
12 CO jZiZq an d J=2-l line ratios for region A2 is more uncertain which may result in 
unreliable radiative transfer results. 

To constrain the a verage temperature and density of each region, we use the radiative 



transfer code RADEX (Ivan der Tak et al.l 120071 ). RADEX assumes a homogeneous medium 



making it a useful tool to constrain the physical conditions such as density and temperature 
using observational data. We created a 3D grid of models varying iV( 12 CO) in steps of 0.1 
x 10 18 cm" 2 from 10 17 to 10 20 cm" 2 , T kin in steps of 1 K from 10 - 1000 K and log(ra(H 2 )) 
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in steps of 0.02 from 2-8. We adopt a 12 CO to 13 CO abundance ratio of 50. We used the 
line widths (dV) from Table 3 to convert N( 12 CO)/dV to N( 12 CO). The RAD EX solutions 
were plotted as T^ in vs. n(H 2 ) with contours for each line ratio ±ler. The solutions are 
picked out where the three line ratios overlap within the ±ler range by inspection by eye 
(Figure HJ). For regions A, A2 and B, a wide range of solutions exist. In particular, there 
are solutions with low (< 30 K) temperature but density only constrained to be > 10 4 cm -3 
and also solutions with lower densities but temperature only constrained to be > 30 K. The 
solutions for regions A, A2, and B suggest that the pressure of the medium is P/k > 10 5 
K cm -3 . Aalto et al. (1997) estimated a temperature of > 50 K and a density of 10 4 -10 5 
cm -3 for region A using 12 CO and 13 CO J=l-0 and HCN J=l-0 line ratios. This is not 
consistent with our RADEX solutions but our analysis consists of more CO lines and we do 
not assume local thermal equilibrium (LTE). The HCN J=l-0 line traces molecular gas at 
a higher critical density than our 12 CO J=3-2 line; howe ver, the HCN lines may be affected 
by IR pumping through a 14 /im vibrational transition ( lAalto et al.l 119951 ) and in extreme 
IR environ ments can be over-luminou s with respect to the lines of other dense molecular 
gas tracers ( jGracia-Carpio et al.l 120061 ) . If the HCN J=l-0 line is over-luminous, radiative 
transfer solutions would show higher densities than the true value. For the sub-regions CI 
and C2, we have better constrained solutions for both temperature and density (Table 5). For 
sub-region C2, the solutions show a slightly wider temperature, density and column density 
range than for region CI. The optical depth of the solutions for sub-regions CI and C2 was 
found to be ~0.5-1.5, ~2.5-4 and ~ 4 for 12 CO J=l-0, J=2-l and J=3-2, respectively. The 
low optical depth (<1) of 12 CO J=l-0 line is likely responsible for the unusual (> 1) 2-1/1-0 
line ratios seen in Arp 299. The assumption of a homogeneous medium adds uncertainty to 
our solutions and the solution column densities may easily be off by a factor of 2 from the 
true values. We stress also that our solutions are averages over the beam. 

Our RADEX solutions for regions A, A2 and B show a wide range of solutions such as 
low (< 30 K) temperature but density > 10 4 cm -3 and also solutions with lower densities 
but temperature only constrained to be > 30 K. One way this can be interpreted is a two 
gas component solution with a moderately dense warm (T]<- m > 30 K) ga s comp onent and a 
dense cold (T^ in ~ 10-30 K) gas component. In Arp 220. iRangwala et al.l ( 1201 ll ) have shown 
that the mid- J to high- J CO transition lines are tracing the warm/hot gas and the low- J CO 
transition lines are tracing the cold gas (T ~ 50 K). The mass of the warm gas is 10% of the 
cold gas mass in Arp 220 (Rangwala et al. 2011). This result suggests that even our 12 CO 
J=3-2 observations may be most sensitive to the dominant colder gas component that our 
RADEX solutions allow. In order to trace the warm/hot ga s component more defin i tively , 
we need to observe Arp 299 in higher 12 CO transition lines. iPereira-Santaella et al.l (120101 ) 
have used mid-infrared Spitzer spectra of IC 694 and NGC 3690 to show that there is 5 x 
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10 8 M & of warm (Tk; n > 100K) H2 gas. Since Arp 299 does have a warm/hot gas component 
up to ~ 500 K, we would expect to see strong emission in higher J-level 12 CO lines. 
The Herschel FTS spectrum for Arp 299 indeed shows a prominent 12 CO transition ladder 
with strong emission up to 12 CO J=13-12 for regions A and B (M. R. P. Schirm; private 
communication) . 



Recently, iPapadopoulos et al.l ( 120121 ) have performed large velocity gr adient (LVG) mod- 
elling of a single gas phase of Arp 299 using single dish CO transition lines. IPapadopoulos et al 



(12012] ) used global 12 CO 



J=3-2 J=4-3 
J=l-0> 7=1-0 



and 



12 co 
13 co 



J=l-0 line ratios and found a best fit of T^ in 



= 30 K and n(Ha) = 10 cm" , well within our RADEX solutions. Other solutions were 
shown to be possible with higher temperatures and low densities, similar to our own RADEX 
solutions for regions A, A2 and B. 



3.2. Arp 299 in Context: Comparison to M33 

LVG modelling has previously been performed on a sample of seven in dividual giant 
molecular clouds (GMCs) in the normal spiral galaxy M33 ( IWilson et al.lll997l ). The average 
of six clouds resulted in a solution of T kin =10-20 K, n(H 2 ) = 4-30 x 10 3 cm" 3 and iV( 12 CO) 
= 3-6 x 10 17 cm -2 . This average solution resembles our RADEX solution for region CI 
except for the beam-a yeraged column density which is significantly lower. Since the sizes of 
the clouds are known ( IWilson fc Scovillelll990[ ) we can determine the average column density 
over the cloud itself using the filling factor of the cloud within the beam, / = -Dcioud/^beam; 
where -D c i oud is the diameter o f the cloud and D h fa m is the physical diameter of the beam 
(85 pc). Using the equation in I Wilson fc Scovilld (119901 ) 



FWHM 



1.2D 



0.5 

pc 



(1) 



where Vfwhm is the line width FWHM of the CO line in km s -1 and D pc is the diameter of 
the cloud in pc and using the average Vfwhm for the six clouds given in Wilson et al. (1997) 
(= 9.3 km s -1 ), we get an average -D c ioud of 60 pc. This gives us a filling factor of / = 0.5 
and a cloud-averaged column density of iV( 12 CO) = 6-12 x 10 17 cm -2 for the M33 clouds. 
Comparing the cloud-averaged column density value to the 12 CO column density of region 
C, we see that the sub-region CI and even C2 have at least 2-3 times the average column 
density of the clouds in M33. The striking thing is that this column density is achieved 
not just in an individual GMC but over a ~1 kpc diameter region. This analysis implies 
that if region C (= CI + C2) is made up of normal GMCs such as those in M33, the beam 
is filled with clouds spaced very close together such that the entire ~1 kpc region has a 
structure similar to a single GMC and with at least 2-3 clouds along each line of sight. This 
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unusual structure may help explain the short gas depletion time and high star formation 
rate observed in region C in particular and Arp 299 as a whole (see § 4). 

We can also calculate the average volume density of region C assuming a spheriodal 
geometry. The average radius in the plane of the sky of region C is R ~ 350 pc. However, a 
diameter of ~ 700 pc along the line of sight would be much greater than the typical disk scale 
height. From the column density comparison, we might expect to encounter a maximum of 
3 GMCs of about 60 pc in diameter. If the space in between the clouds is negligible then 
the thickness of region C is about 200 pc. Using the molecular gas mass (§ 4 and Table 6) 
and < n(H 2 )> = 3MH 2 /47rmH 2 -R 3 we derive an average volume density of 130 cm -3 . We can 
compare this average volume density to the average value of the six clouds in M33. Using 
the equation given in Wilson & Scoville (1990) 

< n (H 2 ) >= 200(D pc /20pc)- 1 (2) 

and the average cloud diameter, we get an average volume density of about 70 cm -3 for the 
clouds in M33. This comparison implies that region C is much denser than a 60 pc cloud in 
M33, even though it covers a much larger spatial scale. 

It has been suggested previously that U/LIRGs have kiloparsec scale medium that is 
like a GMC as a whole (e.g. Downes & Solomon 1998). An alternative possibility that could 
reproduce the observed large column density is that U/LIRGs contain a medium where the 
mean density of individual clouds and clumps is much higher than typical GMCs (perhaps 
similar to the circumnuclear disk in the Galactic center i.e. Christopher et al. 2005). However 
this second situation would result in a radiative transfer solution volume density ('h.radex) 
much greater than that of the clouds in M33 (tilvg(M33)). For region C, we see that 



^radex ~ ^lvg(M33) ruling out the alternative suggestion. Ilmanishi fc Nakanishil (120061 ) 
have found that interferometric HCN/HCO+ J=l-0 line ratios of each region are less than 1 
likely indicating that there is a lack of very dense gas, which also supports our results. This 
is the first direct evidence from a radiative transfer analysis using spatially resolved data of 
a kiloparsec scale structure having the properties of a GMC in a U/LIRG. 



3.3. CO-to-H 2 Conversion Factor in Region C of Arp 299 

Since region CI has a better constrained radiative transfer solution, we can examine 
both our assumed 12 CO abundance (x co = 3 x 10~ 4 ) and also the value of the CO-to-H 2 
conversion factor. We do this by comparing the RADEX solution's 12 CO column density to 
the average H 2 column density within the beam derived from the 12 CO intensity. The CO 
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luminosity of each region can be calculated using the equation 



V 



K km s 1 pc 2 



3.2546 x 10 7 ( 



Jy km s- iyv Mpc' v GHz' 



:i+z) 



(3) 



where Sqo is the flux of the CO line in Jy km s -1 , is the luminosity distance in Mpc, v a 
is the rest frequency of the CO line in GHz and z is the redshift (Paper I). Since the 12 CO 
J=l-0 map does not include short spacings and we do not know how much flux is missing in 
each region, we use the short spacing corrected 12 CO J=2-l map to measure the gas masses 
of each region. We use an average ratio value of 12 CO jZ^Zl =1-4 (Table 4) to estimate the 
gas mass, which is then given by 



M(H 2 ) = 0.57L' ( 



where L' r 



CO(2-l) 



CO(2-l) 



is giv en by Equation [5] and we assume a 



the 12 CO J=l-0 line flDownes k Solomon! Il998h . 



0.8 Mq (K km s" 1 pc 2 ) 



(4) 



for 



Using the short spacing corrected 12 CO J=2-l map at the degraded resolution, the peak 
intensity of region CI is 170 Jy/beam km s _1 or 2.2 x 10 8 K km s _1 pc 2 which corresponds 
to a mass of Mh 2 = 1.3 x 10 8 M using Equation 4 . With a beam size of 3.6" x 2.4" 
this corresponds to an average H 2 column density of 1.7 x 10 22 cm -2 . Using the 12 CO 
column density obtained from the RADEX solution, we get a 12 CO abundance of x co = 
iV( 12 CO)/iV(H2) = 0.6 - 2.4 x 10" 4 . This is 1.3 - 5 tim es lower than the abundance ratio of 
3 x 10~ 4 that is typically assumed (IWard et al.l 120031 ). This difference may be interpreted 
in two ways. One possibility is that the abundance of 12 CO to is in fact lower in this 



region of Arp 299 than the average Galactic value. iRelano et al.l (120071 ) have reported a 
metallicity that is near-solar for Arp 299 which would suggest that the abundance of 12 CO 
to H 2 s hould not differ g reatly from the value that is typically used. I f we use the updated 
work by iDickmanl (119781 ) for "dark clouds" reported in lMagnani et al.l (11988 ). iV(H 2 ) = (4 ± 
2)x 10 5 A( 13 CO) and our assumed [ 12 CO]/[ 13 CO] of 50, we get an abundance of Zi 2 co = (l-3 
± 0.7) x 10 -4 . This value agrees well with our analysis. A second possibility is that the 
12 CO abundance in Arp 299 is the same as in the Galaxy but the value of the CO-to-H 2 
conversion factor, a co , is smaller than 0.8 M (K km s _1 pc 2 ) -1 . If we use the peak intensity 
value of the 12 CO J=2-l map and the mass within the beam derived from the column density 
solution of region CI we get 



"00(1-0) 



3 x 10~ 4 

0.4 ± 0.3( )M (K km s" 1 pc 2 )" 1 . 



(5) 



To get aco to agree with Downes & Solomon (1998), the x co abundance factor must be ~1.5 
x 10~ 4 , lower than our assumed value by a factor of 2. There is no way of telling which 
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scenario is right and it may well be a combination of the two. iPapadopoulos et al.l (120121 ) 
have also found a CO-to-H 2 factor over their global LVG solutions to be 0.35 - 0.42 M (K 
km s _1 pc 2 ) -1 . This range of values agrees well with the value we estimated for region CI 
using the RADEX solutions. Note that our result is independent of the method used by 
Downes & Solomon (1998). 

We can attempt to place a second limit on the CO-to-H 2 conversion factor using the 
dynamical mass of regio n C (=C1+C2). The dy namical mass (Table 3) can be estimated 



using the equation from I Wilson fc Scovilld (119901 ) 



M dyn = 99A^ F 2 WHMJ D(pc)M (6) 

where AVfwhm is the CO line FWHM in km s _1 and -D(pc) is the diameter of the region in 
pc. Note equation [6] is valid only for regions that are gravitationally bound. If we assume 
that the dynamical mass is entirely due to molecular gas, then aco = Afdyn/ L'co- Using the 
total flux for region C from the short spacing corrected 12 CO J=2-l map and a 12 CO 2-1/1-0 
average ratio of 1.4 (Table 4), we find that ctco = 1-4 M & (K km s _1 pc 2 ) -1 . However, the 
dynamical mass is likely greater than the true molecular gas mass as stars may contribute 
some mass and the region as a whole may not be gravitationally bound. Thus, this value of 
«co is strictly speaking an upper limit. 



4. Star Formation Rates and Depletion Times 



We estimate th e far-infrared luminosity of each region using Lfir = 5.5 x 10 L & 
(ISanders et al.l 120031 ) for Arp 299 and the fractions of the total luminosity for each region 
from Alonso-Herrero et al. (2000) based on ground-based mid-infrared observations. From 
the relative fluxes at 10, 12, 20 and 30 /im through a 6" diameter, Alonso-Herrero et al. 
(2000) estimated that region A has about 50%, region B has about 27% and region C has 
about 13% of the total infrared luminosity. The final 10% is believed to be associated with 
H II regions. The star formation rate can be estimated from the infrared luminosity (8-1000 
fim) using the equation from Kenn icutt (1998), M sfr(M^ yr" 1 ) = 4.5 x 10~ 44 LiR(erg s _1 ). 
Kennicutt (1998) used a Salpeter (jSalpeterl Il955l ) initial mass function (IMF) to estimate 
the star formation rate; a double power law IMF with a slope of -1.3 in the range of 0.1-0.5 
Mq and a slope of -2.3 in the range of 0.5-120 M & is more commonly used (i.e. Calzetti et 
al. 2007). The star formation rate is 1.59 times lower with the double power law IMF than 
with the Salpeter IMF. To minimize any contributi on to the inf r ared l uminosity from any 
AGNs present in Arp 299, we adopt the method of iGenzel et al.l (120101 ). who scale the far- 
infrared luminosity by a factor of 1.3 to estimate the total infrared (8-1000 /im) luminosity 
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associated with star formation. This method avoids using the directly measured mid-infrared 
luminosity that is more likely to be contaminated by an AGN. The correction facto r of 1.3 
is an average L IR /L FIR in nearby star-forming galaxies ( IGracia-Carpio et al.ll2008l ). With 
these changes, the star formation rate becomes 



Msfr(M yr 



1.3 x 10 



-10; 



J FIR 



(£©)■ 



(7) 



Table 6 summarizes the derived properties of Arp 299. 



Using the total molecular gas mass and star formation rates of each region, we can 
estimate the amount of time it will take for all the molecular gas to be depleted by star 
formation using 

1.36M(H 2 ) 



^depl 

M SFR 

where the factor 1.36 is to account for the helium gas. The dep letion times of ea ch re- 
gion (Table 6) are very short compared to normal spiral galaxies ( ILeroy et al.l 120081 ) . The 
timescales are more similar to the lifetime of GM Cs i n the Large Magelle nic Cloud (20-30 
Myr; Kawamura et al., 2009). iDaddi et al.l ( 120101 ) and iGenzel et al.l (120101 ) studied the star 
formation laws in disk galaxies and starbursts at low and high redshifts. They provide ev- 
idence that disks and starbursts occupy distinct regions in the molecular gas mass (Mh 2 ) 
versus star formation rate (or Ljr) plane. Starbursts are seen to have 10 times higher Ljr 
at fixed Mh 2 . U/LIRGs have depletion times on t he order of 10 7 yea rs while spiral galaxies 
have depletion times in the range of 0.6 - 2.4 Gyr (IDaddi et al.ll2010l ). 



To increase the gas depletion times in Arp 299 to values similar to normal spiral galaxies, 
either the molecular gas mass is underestimated or the star formation rate is overestimated. 
The molecular gas mass must be increased by 2 orders of magnitude to achieve similar de- 
pletion times to those seen in normal spiral galaxies. However, even if we take the dynamical 
masses to be entirely molecular gas, we do not get close to increasing the gas mass by this 
amount. If the far-infrared luminosity is heavily contaminated by an AGN then the star 
formation rates could be overestimated. If we take the extreme case that an AGN con- 
tri butes 80% to t he tot al infrared luminosity, the largest contribution found for any U/LIRG 
by IGenzel et al.l ( 119981 ) , the depletion timescales in the two nuclei would only be increased 
by a factor of 5. In either case, taking the molecular gas mass or the star formation rates 
to the extreme limits does not give us depletion timescales similar to those of normal spiral 
galaxies. 

The global gas depletion time of Arp 299 suggests that all of the molecular gas will be 
used up in star formation in about half of the merger time (~ 60 Myr; Hibbard & Yun 1999). 
If the molecular gas is not replenished over time, Arp 299 may not reach the ULIRG phase. 
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However, Hibbard & Yun (1999) have detected a tidal tail containing 3.3 x 10 9 M of HI. 
If the gas inflow rate is sufficiently fast then there could be replenishment of molecular gas 
for further star formation. 



5. Conclusions 

We have analyzed the interferometric CO observations of the local LIRG Arp 299 first 
published in Wilson et al. (2008). CO line ratios have been used to constrain the physical 
conditions of each region of Arp 299 using the radiative transfer code RADEX. Our RAD EX 
solutions for regions A, A2 and B show two situations, solutions with low (< 30 K) temper- 
ature but density only constrained to be > 10 4 cm~ 3 and also solutions with lower densities 
but temperature only constrained to be > 30 K. This analysis allows for the possibility of a 
two component gas, but higher J-level 12 CO transition lines are required to trace the warm 
gas conclusively. 

The sub-region CI corresponding to the peak intensity overlap region of Arp 299 shows 
a better-constrained RADEX solution with T^ n ~ 10-50 K, n(H2) ~ 1 - 30 x 10 3 cm" 3 and 
iV( 12 CO) ~ 1-4 x 10 18 cm -2 . This solution is comparable to that of GMCs in the normal 
spiral galaxy M33. The column densities of the clouds of M33 are at least 2-3 times lower 
than the column density for region C (= CI + C2) which may indicate that clouds are spaced 
close together in region C. We estimate the CO-to-H 2 factor, a co to be 0.4 ± 0.3 (3 x 10~ 4 / 
x co ) M (K km s _1 pc 2 ) -1 for region CI, which agrees well with Downes & Solomon (1998) 
for more compact and luminous ULIRGs. Using the dynamical mass of region C (= CI + 
C2), an upper limit to aco was determined to be 1.4 M (K km s _1 pc 2 ) -1 

We have used the CO observations to calculate the gas and dynamical masses of each 
region and infrared data to derive the star formation rate of each region. We find that the 
gas depletion times of each region (20-60 Myr) are about 2 orders of magnitude lower than 
those for normal spiral galaxies and resemble the 30 Myr depletion timescales of GMCs. 
These short depletion timescales may be explained by a higher concentration of dense gas 
than in normal spiral galaxies. If the molecular gas is not replenished over time, Arp 299 
will run out of fuel for star formation in about half the time it will take for the nuclei to 
finally merge. 

The Submillimeter Array is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by 
the Smithsonian Institution and the Academia Sinica. The James Clerk Maxwell Telescope 
is operated by The Joint Astronomy Centre on behalf of the Particle Physics and Astronomy 
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search, and the National Research Council of Canada. We thank the anonymous referee for 
a very useful referee report and S. Aalto for giving us the OVRO CO J=l-0 map. C.D.W. 
acknowledges support by the Natural Science and Engineering Research Council of Canada 
(NSERC). 




J2000 Right Ascension J2000 Right Ascension 

Fig. 1. — Integrated intensity maps of Arp 299. All SMA maps have been corrected for the 
primary beam. Dashed contours indicate negative flux: (top left) 12 CO J=l-0 OVRO map 
from data cube given by Aalto et al. (1997) with contours corresponding to -4, 2, 6, 10, 
20, 30, 40, 50 x 0.9 Jy/beam km s _1 (=ler), (top right) 12 CO J=2-l map with contours 
corresponding to -4, 2, 6, 10, 20, 30, 40, 50 x 1.7 Jy/beam km s" 1 (=la), (bottom left) 13 CO 
J=2-l map with contours corresponding to -2, 2, 4, 6, 8, 10 x 1.4 Jy/beam km s _1 (=lcr), 
(bottom right) 12 CO J=3-2 map with contours corresponding to -4, 2, 6, 10, 20, 30, 40, 50 
x 3.5 Jy/beam km s _1 (=la). Note: The SMA maps from Paper I have been reprocessed 
using CASA. 



REFERENCES 

Aalto, S., Booth, R. S., Black, J. H., & Johansson, L. E. B. 1995, A&A, 300, 369 
Aalto, S., Radford, S. J. E., Scoville, N. Z., & Sargent, A. I. 1997, ApJ, 475, L107+ 



- 15 - 



Alonso-Herrero, A., Rieke, G. H., Rieke, M. J., & Scoville, N. Z. 2000, ApJ, 532, 845 
Armus, L., Heckman, T., & Miley, G. 1987, AJ, 94, 831 
Barnes, J. E. & Hernquist, L. E. 1991, ApJ, 370, L65 
Calzetti, D. et al. 2007, ApJ, 666, 870 

Casoli, F., Willaime, M.-C, Viallefond, F., & Gerin, M. 1999, A&A, 346, 663 

Christopher, M. H., Scoville, N. Z., Stolovy, S. R., & Yun, M. S. 2005, ApJ, 622, 346 

Currie, M. J., Draper, P. W., Berry, D. S., Jenness, T., Cavanagh, B., & Economou, F. 2008, 
in Astronomical Society of the Pacific Conference Series, Vol. 394, Astronomical Data 
Analysis Software and Systems XVII, ed. R. W. Argyle, P. S. Bunclark, & J. R. Lewis, 
650-+ 

Daddi, E. et al. 2010, ApJ, 714, LI 18 

Dickman, R. L. 1978, ApJS, 37, 407 

Downes, D. & Solomon, P. M. 1998, ApJ, 507, 615 

Gehrz, R. D., Sramek, R. A., & Weedman, D. W. 1983, ApJ, 267, 551 

Genzel, R. et al. 1998, ApJ, 498, 579 

-. 2010, MNRAS, 407, 2091 

Gracia-Carpio, J., Garcia-Burillo, S., Planesas, P., & Colina, L. 2006, ApJ, 640, L135 

Gracia-Carpio, J., Garcia-Burillo, S., Planesas, P., Fuente, A., & Usero, A. 2008, A&A, 479, 
703 

Hibbard, J. E. & Yun, M. S. 1999, AJ, 118, 162 

Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 

Imanishi, M. & Nakanishi, K. 2006, PASJ, 58, 813 

Iono, D., Wilson, C. D., Yun, M. S., Baker, A. J., Petitpas, G. R., Peck, A. B., Krips, M., 
Cox, T. J., Matsushita, S., Mihos, J. C, & Pihlstrom, Y. 2009, ApJ, 695, 1537 

Kawamura, A., Mizuno, Y., Minamidani, T., Filipovic, M. D., Staveley-Smith, L., Kim, S., 
Mizuno, N., Onishi, T., Mizuno, A., & Fukui, Y. 2009, ApJS, 184, 1 



- 16 - 



Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 

Leroy, A. K., Walter, F., Brinks, E., Bigiel, F., de Blok, W. J. G., Madore, B., & Thornley, 
M. D. 2008, AJ, 136, 2782 

Magnani, L., Blitz, L., & Wouterloot, J. G. A. 1988, ApJ, 326, 909 

Menendez-Delmestre, K., Blain, A. W., Smail, I., Alexander, D. M., Chapman, S. C, Armus, 
L., Frayer, D., Ivison, R. J., & Teplitz, H. 2009, ApJ, 699, 667 

Murphy, Jr., T. W., Armus, L., Matthews, K., Soifer, B. T., Mazzarella, J. M., Shupe, D. L., 
Strauss, M. A., & Neugebauer, G. 1996, AJ, 111, 1025 

Papadopoulos, P. P., van der Werf, P., Xilouris, E., Isaak, K. G., & Gao, Y. 2012, ArXiv 
e-prints 

Pereira-Santaella, M., Alonso-Herrero, A., Rieke, G. H., Colina, L., Diaz-Santos, T., Smith, 
J.-D. T., Perez-Gonzalez, P. G., & Engelbracht, C. W. 2010, ApJS, 188, 447 

Rangwala, N., Maloney, P. R., Glenn, J., Wilson, C. D., Rykala, A., Isaak, K., Baes, M., 
Bendo, G. J., Boselli, A., Bradford, C. M., Clements, D. L., Cooray, A., Fulton, T., 
Imhof, P., Kamenetzky, J., Madden, S. C, Mentuch, E., Sacchi, N., Sauvage, M., 
Schirm, M. R. P., Smith, M. W. L., Spinoglio, L., & Wolfire, M. 2011, ApJ, 743, 94 

Relano, M., Lisenfeld, U., Perez-Gonzalez, P. G., Vilchez, J. M., & Battaner, E. 2007, ApJ, 
667, L141 

Salpeter, E. E. 1955, ApJ, 121, 161 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A., & Soifer, B. T. 2003, AJ, 126, 
1607 

Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 

Sanders, D. B., Scoville, N. Z., Young, J. S., Soifer, B. T., Schloerb, F. P., Rice, W. L., & 
Danielson, G. E. 1986, ApJ, 305, L45 

Sanders, D. B., Soifer, B. T., Elias, J. H., Madore, B. F., Matthews, K., Neugebauer, G., & 
Scoville, N. Z. 1988, ApJ, 325, 74 

Sargent, A. & Scoville, N. 1991, ApJ, 366, LI 

Schirm, M. R. P. 2011, private communication 



-17- 

van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, 
A&A, 468, 627 

Ward, J. S., Zmuidzinas, J., Harris, A. I., & Isaak, K. G. 2003, ApJ, 587, 171 
Wilson, C. D. & Scoville, N. 1990, ApJ, 363, 435 
Wilson, C. D., Walker, C. E., & Thornley, M. D. 1997, ApJ, 483, 210 
Wilson, C. D. et al. 2008, ApJS, 178, 189 



This preprint was prepared with the A AS IATgX macros v5.2. 



-18 - 




H^B^SS 8 34 s 33° 32 s 31' 

J2000 Right Ascension 



Fig. 2. — Feathered (short spacing corrected) SMA maps: {top) 12 C0 J=2-l map with 
contour levels as in Figure [TJ The light blue lines indicate the area where the flux was 
integrated for each region (see Table 6); (bottom) 12 CO J=3-2 map with contour levels as in 
Figure [TJ 
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J2000 Right Ascension J2000 Right Ascension 

Fig. 3.— Line ratio maps: {Top Left) ^§g J=2-l; {Top Right): 12 CO ; {Bottom Left): 

SMA only 12 CO j^fzf ; {Bottom right): Feathered 12 CO j=§Ef ratio map. Only emission 
that is > 2a in both maps is included. 
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log(n(H 2 )) (crrT 3 ) log(n(H 2 )) (cm J ) 

Fig. 4.— Sample RADEX plots of each region of Arp 299: (top left) Region A for iV( 12 CO) 
= 2 x 10 18 cm - 2 , (top middle) Sub-region C2 for iV( 12 CO) = 3 x 10 18 cm " 2 , (top left) 
Sub-region CI for iV( 12 CO) = 2 x 10 18 cm ~ 2 , (bottom left) Region A2 for iV( 12 CO) = 4 
x 10 18 cm ~ 2 and (bottom right) Region B for N( 12 CO) = 2 x 10 18 cm ~ 2 . The thick solid 
lines denote solutions. The entire range of solutions is not shown. 



Table 1. Interferometric Data For Arp 299 



CO Transition Line Total Flux a Beam rms 





(Jy km s x ) 




CO 




(mJy/beam) b 


12 C0 J=1-0 C 


377± 9 


2.5 


X 


2.2 


9 


12 C0 J=2-l 


1870 ± 30 


3.0 


X 


1.8 


17 


13 C0 J=2-l 


40 ± 2 


3.2 


X 


1.8 


10 


12 C0 J=3-2 


2820 ± 40 


2.2 


X 


1.9 


33 



a measurement uncertainty only; calibration uncertainty is 
20% (Paper I). 

b rms values for a 20 km s _1 channel width for all lines except 
13 CO J=2-l which uses 50 km s" 1 channel width. 



c Map obtained from lAalto et al 
Valley Radio Observatory. 



(119971 ) ; observed with Owens 
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Tablc 2. Comparision of JCMT and Combined (Feathered) Maps a 



Total Flux (Jy km s" 1 ) 13 



Transition SMA Feather JCMT Missing c Recovered d 

( 12 CO) (%) (%) 



J=2-l 1870 ± 30 2380 ± 30 2320 ± 120 20 ± 30 100 ± 30 
J=3-2 2820 ± 40 6350 ± 50 8770 ± 120 70 ± 30 70 ± 40 



al3 CO J=2-l and 12 CO J=l-0 single dish maps are not available 

b measurement uncertainty only; calibration uncertainty is 20% for SMA 
and the JCMT 12 CO J=2-l maps and 30% for the JCMT 12 CO J=3-2 (Paper 

I)- 

c Missing flux in the SMA only map relative to the JCMT map. Uncertainty 
calculation includes calibration uncertainty. 

d Recovered flux in the feathered map relative to the JCMT map. Uncer- 
tainty calculation includes calibration uncertainty. 



Table 3. 



Region 


Velocity FWHM a 


Deconvolved Source Diameter* 3 


M dyn 




(km s _1 ) 






(PC) 


(10 9 M ) 


A (IC 694) 


325 


(1.7 ± 0.5) x 


(1.2 ± 1.0) 


380 x 270 


(3.4 ± 2.1) 


A2 (IC 694 disk) 


90 










B (NGC 3690) 


185 


(2.0 ± 0.6) x 


(1.4 ± 1.0) 


450 x 310 


(1.3 ± 0.7) 


CI (Overlap) 


80 


(2.5 ± 0.8) x 


(1.6 ± 1.3) 


560 x 360 


(0.3 ± 0.2) 


C2 (Overlap) 


80 


(4.3 ± 2.1) x 


(1.5 ± 1.4) 


960 x 330 


(0.4 ± 0.3) 



a Uncertainty in velocity FWHM is 10 km s 1 
b Beam size is 2.2" x 1.9" 

c Region A2 cannot be fit with a Gaussian because it is too close to region A 
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Table 4. CO Line Ratios of Arp 299 a 



Region 


12 CO f 3 ~ 2b 


-co ^§Ef 


i2nr\ J=2-i 

UU ,7=1-0 


12co J-2-1 
13 co J— z 1 


A 


0.74 ± 0.21 


0.72 ± 0.20 


1.37 ± 0.39 


30.9 ± 4.4 


A2 


0.76 ± 0.22 


0.58 ± 0.16 


1.32 ± 0.38 


12.6 ± 1.8 


B 


0.73 ± 0.21 


0.73 ± 0.21 


1.66 ± 0.47 


27.3 ± 3.9 


CI 


0.52 ± 0.15 


0.51 ± 0.15 


1.24 ± 0.35 


17.2 ± 2.4 


C2 


0.65 ± 0.19 


0.60 ± 0.17 


1.23 ± 0.35 


11.6 ± 1.6 



a All ratios are in terms of integrated brightness temperatures 
at the peak intensity value found in the 12 CO J=3-2 map. Uncer- 
tainties include calibration errors only. All maps were degraded 
to a resolution of 3.6" x 2.4". 

b Using short spacing corrected maps 

c Using SMA only maps 



Table 5. Radex Results 



Region 


N( 12 CO) 




n(H 2 ) 




(cm -2 ) 


(K) 


(cm -3 ) 


A 


1-8 x 10 18 


10-500 


>io 2 - 5 


A2 


2-30 x 10 18 


20-500 


>io 2 - 5 


B 


1-6 x 10 18 


10-1000 


>io 2 - 5 


CI 


1-4 x 10 18 


10-50 


10 3 -10 45 


C2 


2-9 x 10 18 


10-200 


10 3 -10 5 
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Table 6. Gas Mass and Depletion Times 



Region 


MSFR 

(M e yr- 1 ) 


12 C0 J=2-l Flux a 
(Jy km s _1 ) 


r be 
-^CO(2-l) 

(10 9 K km s" 1 pc 2 ) 


M(H 2 ) C 
(10 8 M ) 


tdepl 

(Myr) 


A 


35.8 


920 ± 100 


1.2 ± 0.2 


7 ± 1 


27 ± 5 


A2 


<7.2 d 


320 ± 20 


0.4 ± 0.1 


2.3 ± 0.5 


>43 


B 


19.3 


390 ± 40 


0.5 ± 0.1 


3.0 ± 0.6 


21 ± 4 


C = CI + C2 


9.3 


510 ± 30 


0.7 ± 0.1 


4.0 ± 0.7 


60 ± 10 


Total 


71.5 


2380 ± 20 


3.1 ± 0.6 


18 ± 4 


34 ± 7 



a measurement uncertainty only; calibration uncertainty is 20% 
b Using D L = 46 Mpc 

Uncertainty is calculated using the calibration uncertainty 

d M SFR calculated assuming the region contains 10% the total infrared luminosity of Arp 
299. 



